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PSC occur, the PGS power-voltage characteristic curve exhibits multiple peak values; that is, the curve 
comprises a global maximum power point and multiple local maximum power points. Current literature 
includes various studies of global maximum power point tracking (GMPPT) algorithms and hardware 
architectures suitable for PSC; because the substantial quantity of PSC literature, this subject must be 
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1. Introduction 


According to a technical report by the European Photovoltaic 
Industry Association (EPIA), more than 31 GW of photovoltaic (PV) 
capacity was installed globally in 2012; aggressive projections 
indicate that this number could increase to 84 GW in 2017 [1]. 
Weather and environmental factors affect the output character- 
istics of PV cells, that is, cell output voltage and current vary 
according to changes in irradiance and temperature. Therefore, a 
unique maximum power point (MPP) exhibits under specific 
irradiance and temperature. The maximum power point tracking 
(MPPT) algorithm is crucial in attaining the maximal PV power, 
facilitating optimal PV cell performance. Researchers have pro- 
posed numerous MPPT algorithms that demonstrate excellent 
tracking efficiency in uniform insolation conditions. However, PV 
generation systems (PGSs) in urban environments installed on 
rooftops or eaves are easily shaded by neighboring buildings, 
clouds, or dust coverage, causing PGSs to generate power output 
that differs from that in ideal insolation conditions. Shading is 
especially problematic regarding PGSs that exhibit wide-range 
layouts primarily because complete shading of a large area is 
unlikely, causing partially shaded conditions (PSC). When PSC 
occur, certain PV cells are exposed to normal, non-shaded condi- 
tions, and output normal power levels; however, other cells yield 
decreased output levels because of shading, generating a mis- 
match in the overall PGS. Because of the characteristics of PV cells, 
shaded conditions significantly reduce cell output current but not 
output voltage. This mismatch in power generation influences 
series-connected PGSs; however, traditional MPPT algorithms are 
designed to address only uniform insolation and cannot ade- 
quately account for PSC. Consequently, developing new MPPT 
techniques suitable for application in PSC is critical. Numerous 
studies regarding this subject have been published; thus, compos- 
ing an accurate and comprehensive reference document that 
details information related to various MPPT techniques applicable 
to PSC is essential. 

Various reviews [2-10] of MPPT methods for PGSs exist in the 
current literature; of these studies, most of the reviews [2-7] 
researched only traditional MPPT methods, and Ishaque and Salam 
[8] and Salam et al. [9] focused on both traditional methods and 
those suitable for PSC. When PSC occur, the power-voltage (P-V) 
characteristic curve of the PV modules becomes complex, exhibit- 
ing multiple peak values. Therefore, traditional MPPT methods are 
not suitable for use in PSC. To provide PGS researchers and 
engineers with useful reference information, this study explored 
global maximum power point tracking (GMPPT) methods suitable 
for application in PSC. First, the basic operating principles of 
various GMPPT methods, including their advantages and limita- 
tions, were described. Regarding firmware-related techniques, 
various indices were evaluated, including whether the design 
parameters were relevant to the system, GMPPT time, GMPPT 
accuracy, algorithm complexity and required sensors. In terms of 
hardware-related techniques, the advantages and disadvantages of 
various architectures were discussed. Finally, an in-depth simula- 
tion on some highly cited GMPPT methods was made and conclu- 
sion was offered accordingly. It was difficult to propose a fair 
benchmark for various techniques mentioned in the literature 
primarily because researchers used distinct PGSs (including power 
rating and power converter topology) and shading patterns to test 
their proposed algorithms. Nonetheless, the authors of the current 
study endeavored to employ several fair evaluation indices and 
verify the correctness of this performance evaluation using the 
simulated results in Section 4. The information gathered herein 
should provide researchers with comprehensive reference infor- 
mation and recommendations for future studies addressing PSC. 
Although Bidram et al. [10] also compared GMPPT methods, this 


study additionally classified various GMPPT methods into cate- 
gories and explored their usage limitations and design considera- 
tions. Furthermore, shading detection methods proposed in the 
literature were compared, and a simulation platform was provided 
for verifying the performance of some highly cited GMPPT meth- 
ods. To the best of the authors’ knowledge, no existing literature 
has provided this type of categorizing and testing. 


2. Exploration of the effects of partially shaded conditions 


When PV cells are connected in series, a bypass diode is 
typically added in parallel for protection, preventing damage due 
to overheating when one or more PV cells are shaded or shorted. 
However, adding a bypass diode changes the current-voltage (J-V) 
characteristic curve when the PGS is shaded by environmental 
conditions such as clouds, tree shade, and building shadows. 
Furthermore, the P-V characteristic curve also changes accordingly 
and exhibits multiple peak values. Current studies on the effects of 
PSC can be divided into three categories: (1) deriving analytical 
mathematical formulas to describe the I-V and P-V characteristic 
curve [11,12]; (2) establishing a simulation or experimental plat- 
form that can evaluate the effects of PSC [13-20]; and (3) evaluat- 
ing the power losses caused by PSC [21-25]. This section explores 
the effects of PSC in PGSs. 


2.1. Basic characteristics of PV cells 


Based on their electrical characteristics, PV cells are equivalent 
to a one-diode model as shown in Fig. 1. This model includes a 
current source Is, a diode current lą, an equivalent parallel 
resistance Rp, and an equivalent series resistance R,. The output 
current of the PV module Ipy is the difference between the 
photogenerated current I; and diode current Ig, as shown in (1): 


_ q(Vpy +IpyRs) Vpv +Rslpv 
Iv = Is—Iso(exp( akT 1 Rp 


where Iso is the saturation current, Vpy is the output voltage of PV 
module, q is the charge of a single electron, a is the ideal factor of 
the diode, k is the Boltzmann constant, T is the temperature in 
Kelvin. 


(1) 


2.2. Effect of partial shade conditions on photovoltaic generation 
systems 


PGSs are formed by several PV modules, which comprise 
several PV cells connected in series or parallel. Parallel connections 
can increase the output current of the module, and serial connec- 
tions can improve the output voltage. When PSC occur, the shaded 
PV cell becomes a load instead of a generator, producing a hot spot 
that can subsequently destroy PV cells. Therefore, commercially 
available PV modules often comprise a parallel-connected bypass 
diode to prevent this problem. When PSC occur, currents flow 
through the bypass diode instead of the shaded PV module, 
generating multiple steps in the I-V characteristic curve and 


Fig. 1. Equivalent circuit of photovoltaic cells. 
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Fig. 2. A 5s1p PGS with various shading patterns. 
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Fig. 3. The characteristic curves of the 5s1p PGS in Fig. 2: (a) P-V characteristic curve; (b) I-V characteristic curve. 


multiple peak values in the P-V characteristic curve. For example, 
Fig. 2 shows a PGS comprising five PV modules connected in series 
(denoted as 5s1p system). Its P-V and I-V characteristic curves for 
different shading patterns are illustrated in Fig. 3(a) and (b), 
respectively. Fig. 3 shows that GMPP under PSC may occur on a 
larger voltage range compared with that in PGSs under uniform 
insolation. 


3. Global maximum power point tracking techniques for 
partial shading conditions 


In this section, GMPPT techniques suitable for PSC are explored. 
The solutions for PSC proposed in current literature can be divided 
into firmware-based [26-69], and hardware architecture-based 
solutions [70-113]. Fig. 4 displays common PGS architectures, 
illustrating that when PGSs use a centralized architecture, only 
one set of power converter is required; therefore, firmware must 
be used to address the multiple peak values that the PSC cause in 
P-V characteristic curves. On the other hand, when PGSs use 
distributed architectures, multiple sets of power converter exist; 
therefore, each power converter can track its own MPPT. Because 


hardware-based solutions relate to the power converter topo- 
logy and system design, this study mainly focuses on exploring 
firmware-based GMPPT techniques. Fig. 5 shows the P-V charac- 
teristic curve of a typical 5s1p system under a specific shading 
pattern; this system was used as the basis to clarify operating 
concepts and design guidelines regarding various GMPPT algo- 
rithms. (Fig. 6). 


3.1. Firmware-based GMPPT techniques 


3.1.1. Two-stage method 

Fig. 5 shows that when PSC occur, the P-V characteristic curve 
of PGS exhibits multiple peak values, and the number of peak 
values correlates with the shading pattern. By dividing this P-V 
characteristic curve into multiple single-peak P-V characteristic 
curves, traditional MPPT methods (e.g., perturb & observe (P&O) or 
incremental conductance (IncCon) methods) can be used to 
determine the maximal value of each single-peak curve. Therefore, 
the two-stage technique involves using methods proposed in the 
first stage to determine the approximate GMPP location (i.e., an 
interval), subsequently using traditional MPPT methods to deter- 
mine the precise GMPP location. Thus, the key to the success of the 
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Fig. 5. P-V characteristic curve under a specific shading pattern. 


two-stage method is the accuracy of the GMPP interval obtained in 
the first stage. The following sections clarify the two-stage method 


used in the literature. 


3.1.1.1. System characteristic curve method. Works in [26-28] 
employed a similar technique, using a preset linear function to 
move the operating point (OP) near the GMPP under PSC, as 
shown in Fig. 7. This linear function is related to various system 
parameters, including open circuit voltage and short circuit 
current. The linear function used in [26] can be expressed as 
follows: 

Vori 
lpm 


Rm= (2) 


Fig. 6. Concept of dividing multiple-peak P-V characteristic curve into several 
single-peak curve. 


The linear function used in [27] can be expressed as follows: 


Pype(t) = [TOIME] (3) 
And the linear function used in [28] can be expressed as follows: 


( rE) x Ipin] (4) 


) x Ip[n ~ ¢ 
out 

The advantage of this method is its tracking speed. Because the 

first stage only requires one step to move the OP to the vicinity of 

the GMPP, the tracking time of this method is similar to that of 

traditional MPPT methods. However, disadvantages exist. First, 


V= 
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Fig. 7. System characteristic curve method. 
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Fig. 8. Incorrect tracking of GMPP using the system characteristic curve method. 
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Fig. 10. Incorrect tracking of GMPP using the method proposed by Patel. 


obtaining the open circuit voltage and short circuit current of the 
system requires open or short circuits. Implementing this in an 
actual system can cause power losses or safety concerns. Second, 
although the linear function used in this method can move the OP 
near the GMPP in certain shading patterns, an erroneous location 
of the new OP can prevent tracking the GMPP; thus, this method 
cannot guarantee that the GMPP can be tracked. This concept is 


shown in Fig. 8. Furthermore, when shading patterns become 
complex, the probability that this method can successfully track 
GMPP decreases. Ref. [29] is one conference paper that used a 
concept similar to this method. 


3.1.1.2. Two-stage searching method. In a two-stage searching 
method, the first stage is used to search for the GMPP location 
interval; traditional MPPT methods are then used in the second 
stage to find the precise GMPP location. By using this method, 
Patel and Agarwal [30] employed various shading pattern 
observations as clues to determine the GMPP location. As shown 
in Fig. 9, the basic search rules for [30] are as follows: 


1) 0.85 Vocan (Voc,an represents the total open circuit voltage of the 
system) is used as the P&O search starting point to find and 
record a peak value. 

2) The peak value found in Step 1 is used as the basis for moving 

the OP one large step (the reference recommended 0.6-0.7 Voc, 

one (Vocone represents the open circuit voltage of a single 
module)) to the left, then conduct P&O to search and record 
the next peak value. 

If the obtained peak value is higher than the previous peak 

value, Step 2 is repeated; if the determined peak value is 

smaller than the previous peak value, then the previous peak 
value is the GMPP. 


Ww 
© 


Ref. [30] is the paper with the highest number of citations 
among similar literature; its method is easy to implement, and can 
be integrated into traditional PGS firmware; however, disadvan- 
tages exist. First, the authors of [30] found that when P-V curve is 
traversed from either side, the magnitude of the peaks increases; 
after reaching the GMPP, the magnitude of the subsequent peaks 
continuously decreases. Therefore, Step 3 can be used to deter- 
mine the GMPP. However, the current study explored all possible 
shading patterns and found that there exists some exception about 
this claim (for detailed descriptions, please refer to Section 4). 
Consequently, this method cannot ensure successful GMPPT. 
Fig. 10 shows one example of this exception. In addition, the 
parameter used in this method (i.e. value of the large step) affects 
the probability of tracking the GMPP. Another disadvantage of this 
method is the tracking speed. Because each LMPP must be 
determined using P&O method, this method requires more track- 
ing steps compared with other methods. 

Works in [31-33] used similar methods, employing a large 
interval to search the entire P-V characteristic curve (called 
global-biased search) to determine the largest peak value. Next, 
a refined search is conducted around the largest peak value found 
in previous stage to locate the GMPP (called local-biased search). 
Of these three approaches, Renaudineau et al. [31] used aq, 
Koutroulis et al. [32] used fixed power intervals, and Kouchaki 
et al. [33] used 0.8 Vocone as the step for global-biased search. The 
concept of this kind of method (Ref. [32] by Koutroulis) is shown in 
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Fig. 11. Method proposed by Koutroulis. 
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Fig. 12. Incorrect tracking of GMPP using the method proposed by Koutroulis. 


Fig. 11. The advantage of this kind of method is its simplicity; it can 
easily be integrated with traditional PGS firmware (requiring only 
an additional global-biased search stage). The disadvantage of this 
kind of method is that the selected global-biased search step 
affects the probability of finding the GMPP. When the selected step 
is large, the system requires comparatively less time to locate the 
GMPP, but may miss the GMPP, as shown in Fig. 12. Conversely, 
when a small step is selected, the system has an increased 
probability of finding the GMPP, but requires relatively more time. 
To use this method, a tradeoff between the search step size and 
the time required must be addressed (for details, please refer to 
Section 4). In addition, the method in [32] used power as the step 
variable for the search. Although this thinking is intuitive and 
yields a comparatively fast tracking speed, this method requires 
using a specifically designed power converter that can accept 
power command for constant power control. Kashif et al. [34] and 
Kouchaki et al. [35] are conference papers used similar two-stage 
searching concepts. 


3.1.2. Segmental search method 

The basic principle of the segmental search method is similar to 
that of the two-stage searching method described in 3.1.1.2. The 
primary difference is that two-stage searching method uses 
traditional MPPT methods (such as P&O) as the basis for the 
second stage search. Segmental search, however, is based on 
mathematical theories, finding the GMPP by gradually reducing 
the search range. 


3.1.2.1. DIRECT method. The theoretical basis for the method used 
in [36] is that the P-V characteristic curve of PV modules conforms 
to Lipschitz characteristics; therefore, the dividing rectangle 
(DIRECT) technique can be used to find the peak values. DIRECT 
refers to dividing the search area into three areas of equal intervals 
by using mathematical methods, subsequently using mathematical 
equations to find potentially optimal intervals (POIs) that exhibit 


potential oe values. The rules for POI search are as follows: 


a a bi) 


fapt k= b) fe p+ vi (5) 


(Gi 


f(q)+K 4+ DD 5 F nant €lf man (6) 


The search area is divided into three equal portions to locate 
the POIs before further subdivision, as shown in Fig. 13. 

Ref. [36] has been frequently cited and compared. The advan- 
tage of this method is that it is based on a solid mathematical 
foundation and facilitates rapid tracking; its disadvantages are 
similar to those of the two-stage searching method. In certain 
shading patterns, this method cannot track the GMPP, although 
this is unlikely. Another potential disadvantage is that although 
the mathematical equations used in this method are relatively 
simple and the segmentation concept can be easily implemented, 
the author added numerous exceptions during actual implemen- 
tation to avoid missing the GMPP. For example, the author 
proposed conducting nine interval searches within the search 
range; however, if the POIs that contained the two largest values 
were not adjacent to each other, then a wide area search with 27 
intervals had to be re-conducted. Conversely, if the obtained POIs 
were adjacent to each other, then segmentation could be con- 
ducted on the POI with the maximal value. This type of exception 
adds complexity and difficulty when firmware engineers imple- 
ment the algorithm based on the DIRECT method. In addition, 
this method cannot be directly integrated into traditional PGS 
firmware. 


3.1.2.2. Fibonacci methods. Works in [37] and [38] used Fibonacci 
sequences as a basis for GMPP searches. Compared with the 
DIRECT method in 3.1.2.1, the primary difference of this method 
is using the Fibonacci sequence as the mathematical basis for 
segmentation, as shown in Fig. 14. The equation is as follows: 


F(n—1) 
X2 =X3 + | FOD (X4—X3) (7) 
A 
im Re 
v 1 
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Fig. 14. The concept of the Fibonacci method. 
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Fig. 15. ANN with three layers (one hidden layer). 
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x1 =x+| las x3) (8) 


where F(n—1) and F(n) represent the (n— 1)th and nth number 
in the Fibonacci sequence. The advantages and disadvantages of 
this method are the same as those in the DIRECT method. 

Ref. [39] and Ref. [40] are conference papers related to the 
segmental search method. Fan et al. [39] used a random number 
method to randomly sample power values under six distinct 
voltages, subsequently using the point with the highest power 
value (of the six selected points) to determine the new search 
range. Repeat random sampling was conducted on the new range 
until convergence was achieved. Escobar et al. [40] used the 
geometric characteristics of the I-V characteristic curve to conduct 
a sweep search with variable steps. 


3.1.3. Soft computing method 

Soft computing (SC) refers to using computers to simulate the 
biochemical processes of natural intelligence systems (e.g., human 
perception, brain structure, and evolutionary process) to solve 
non-deterministic polynomial-time hard (NP-hard) or NP- 
complete problems. Because finding the GMPP in multi-peak P-V 
characteristic curves caused by PSC can be viewed as an optimiza- 
tion problem, scholars have used SC techniques to address such 
problems. SC methods include artificial neural networks (ANNs), 
fuzzy logic control (FLC), evolutionary computation (including 
genetic algorithms (GAs), differential evolution (DE), particle 
swarm optimization (PSO), ant colony systems (ACSs)), and chaotic 
search (CS). 


3.1.3.1. Artificial neural network. An ANN is a computational model 
that uses interconnected artificial neurons to imitate the ability of 
biological neural networks. An ANN with three layers is shown in 
Fig. 15. From Fig. 15, users can input information (such as 
irradiance level and the temperature of PV modules) into the 
input layer and calculate the output (such as the GMPP location) 
according to the weighting values of ANN. For example, the input 
and output relation of an ANN as shown in Fig. 15 can be obtained 
using (9) to (12). The neurons in the input layer are responsible 
only for transmitting the input information to the hidden layer. 
The input of the number j unit in the hidden layer can be 


expressed as follows: 


2 
ni = > will; +b; (9) 
i= 


where wj is the weighting of the ith input unit connected to the 
jth hidden layer unit, and b} represents the bias value of the jth 
unit in the hidden layer. The output of the jth unit in the hidden 
layer can be expressed as follows: 


Vj =f, (n}) (10) 


where f; is the activation function used in this layer. Commonly 
used activation functions include liner function, hyperbolic tan- 
gent sigmoid, and Log-sigmoid. Notably, if a low-cost microcon- 
troller with low calculation capability is used to implement the 
ANN, the linear function is preferred because it is relatively easy 
to apply. 

The input of the k™ unit in the output layer can be expressed as 
follows: 


2 10 2 
m= DY Wyuj+ by (11) 
im1 


where w is the weighting value of the jth hidden layer unit 
connected to the kth output unit, and bz represents the bias value 
of the kth output unit. The output of the kth output unit can be 
expressed as follows: 


Vk = f(n?) (12) 


Similarly, fə is the activation function used for the output layer. 

By using large amounts of training data, the ANN continually 
adjusts the weighting and bias values, allowing the network- 
calculated output to approximate the target output. One of the 
most common ANN learning methods is the back-propagation 
method. The advantage of ANNs is their parallel computing 
capability. Other SC methods may require multiple iterations to 
obtain the optimal solution, whereas ANNs can use simple multi- 
plication and addition to rapidly calculate output. Therefore, ANNs 
enable rapid calculation. However, the accuracy of an ANN is 
determined based on its training data. If the training data are 
insufficient, or the data do not cover the entire problem space, 
then the accuracy of the ANN will be reduced accordingly. 
Currently, the ANN-based MPPT method can only be used in 
uniform insolation conditions [41,42] primarily because at a single 


Y.-H. Liu et al. / Renewable and Sustainable Energy Reviews 41 (2015) 436-453 443 


irradiance level, the location of the MPP is only related to the 
irradiance and the temperature. Thus, the training data are simple. 
When PSC occur, the irradiance of the module, module tempera- 
ture and shading pattern all affect the MPP location. Consequently, 
the training data required by the ANN substantially increase, and 
these data are not easily collected. This is the primary limitation of 
applying ANNs to GMPP searches. Therefore, only works in [43] 
and [44] in current literature proposed using ANNs as the primary 
method for addressing PSC. The input variables used in [43] and 
[44] were the averaged irradiance of selected modules, and the 
design was affected by the arrangement method of PV modules. 
Therefore, if the architecture of the PV array changes (e.g., adding 
new PV panels), the ANN must be re-trained. In addition, irra- 
diance and temperature sensors are more expensive compared 
with the voltage and current sensors used in other GMPP methods. 
Ref. [45] is one conference paper that used ANN concepts. 


3.1.3.2. Fuzzy logic control. Traditional control system design 
requires understanding the system being controlled; that is, using 
precise mathematical models to describe the system. However, when 
the controlled system becomes overly complex, it is often difficult to 
use system identification method to establish a system model. On the 
other hand, FLC uses fuzzy set theory to convert the linguistic values 
directly into automatic control strategies using expert knowledge or 
operator experience. The advantages of FLC include the following: 1) 
no precise mathematical model is required; and 2) FLC can integrate 
the knowledge of human experts into the controller design process. 
When conducting fuzzy control, the behavior of the system being 
controlled is described by a set of fuzzy rules. These fuzzy rules use 
semantic fuzzy information rather than mathematical equations. 
Therefore, FLC can convert the knowledge of human experts into 
fuzzy control rules, reducing the complexity of control system design. 
Fig. 16 shows a typical block diagram of FLC. From Fig. 16, a 
FLC comprises three primary parts: fuzzification stage, a fuzzy 
inference engine, and defuzzification stage. In the fuzzification 
stage, membership function is used to convert the input data into 
linguistic values. The fuzzy inference engine must be used in 
conjunction with the fuzzy rule base, and can be used to calculate 
the required control quantity. In the defuzzification stage, 
mathematical equations are used to convert the output variable 
from a linguistic value to a crisp value. The center of gravity and 
center of area methods are often used for defuzzification. Regarding 
FLC, establishing a rule base is critical in the design. The if-then rules 
formed using linguistic variables can be employed to establish the 
FLC system input and output relations, and must be designed to 
agree with the knowledge of experts concerning a specific problem. 
FLC is particularly suitable for use in non-linear, time-varying 
systems and systems with incomplete models. Thus, FLC can be 
used to deal with PSC problems. Karatepe et al. [46] used distributed 
architecture and utilized FLC to replace traditional MPPT methods. 
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Fig. 16. Basic FLC framework. 


Because each power converter contains a designated MPPT controller, 
the system guarantees that the GMPP can be tracked. The advantage 
of the method used in [46] is its rapid tracking speed and high 
tracking accuracy; its disadvantage is that it exhibits a higher 
hardware cost compared with that of centralized architectures. 
Alajmi et al. [47] used a method similar to the two-stage searching 
method in 3.1.1.2 to conduct GMPPT. This method first swept the P-V 
characteristic curve and recorded various LMPPs, then replaced P&O 
by using FLC to conduct GMPPT. This method yields a rapid tracking 
speed, but demonstrates disadvantages similar to those of other two- 
stage searching methods. Ref. [48] is one conference paper that used 
FLC concepts. 


3.1.3.3. Genetic algorithm. A GA is an evolutionary computation 
(EC) method that mimics the “survival of the fittest” evolutionary 
law in the natural world; its primary operators are reproduction, 
crossover, and mutation. The process of using GA to obtain the 
optimal solution for a problem is as follows: 1) the search 
parameters are coded as binary strings in a chromosome; 2) this 
process is randomly repeated to produce N initial species (strings); 
3) the fitness function is designed based on the solution 
conditions. Species that exhibit high fitness values are selected 
for the mating pool (i.e., the reproduction process); and 4) the 
crossover and mutation processes are computed to complete one 
generation of GA. This process is repeated to produce the species 
with the highest fitness value, as shown in the evolutionary 
process in Fig. 17. GA has been used to optimize the parameters 
of other algorithms (such as FLC or ANNs); only Shaiek et al. [49] 
used GA to solve PSC problems. Compared with other EC methods, 
implementing GA is complex and difficult to achieve using low- 
cost microcontrollers. Consequently, Ref. [49] only provided 
simulated results. Ref. [50] is one conference paper that used GA 
concepts. 


3.1.3.4, Differential evolution. DE is also an EC method, and uses 
architecture similar to that of GA. The primary difference is that DE 
uses special differential operators to replace the crossover process 
in GA and produce the next generation. Comparing to GA, DE is 
simple to implement and requiring few parameter adjustments. 
Currently, no journal paper has used DE to solve PSC problems. 
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Fig. 17. The evolutionary process of GA. 
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Ref. [51] and Ref. [52] are conference papers based on DE, and they 
provided only simulated results. 


3.1.3.5. Ant colony system. An ACS is an EC method based on swarm 
intelligence; its primary advantage is immediately adapting command 
values according to environmental changes. Thus, ACSs are suitable for 
conducting MPP tracking in changing environments. In an ACS, each 
agent selects its path randomly at first. If the path the agent chooses is 
short (has high fitness value), the agent leaves concentrated 
pheromone on the path. In the next iteration, the agent chooses its 
path based on the concentration of pheromone on that path. The more 
concentrated the pheromone is, the higher probability is that the 
agent chooses that path. The probability of the agent moving from 
position i to position j is as follows: 


p 
k Tjj O 13 
i a aR a3) 
UO 
s e allowed, 


where 7; is the pheromone quantity on path (ij), nij is the 
initial value of pheromone on path (ij), & and B are adjustable 
parameters, and the updated pheromone equation can be 
expressed as follows: 


At = 5 Ati (14) 


k=1 


T(t + n)= pti(t) + At; 


In (14), Ark represents the quantity of pheromone that the k" 


agent left on path (ij) in this iteration, whereas Ti represents the 
total pheromone quantity left on path (ij). p is a coefficient such 
that (1- p) represents the evaporation rate. The frequently used 
Ati calculation method is as follows: 


k E 
At; = 4 
ij k 


In (15), Q is a constant and Lẹ represents the distance of the 
path that the k agent passes in this iteration. 

Based on the principle represented by the mentioned equa- 
tions, Jianga et al. [53] used ACS to solve GMPPT problems. 
However, ACSs are difficult to implement using low-cost micro- 
controllers; therefore, Ref. [53] provided only simulated results. 


if the k—th ant passes path (i,j) in this iteration 
, other 


(15) 


3.1.3.6. Particle swarm optimization. PSO is also an EC method 
based on swarm intelligence. In the PSO method, each particle is 
defined by its own position and velocity. The behavior of particles 
within the swarm is influenced by the experiences of neighboring 
particles. Each particle follows the current best-performing particle to 
search within the solution space. In the PSO method, particles are 
initially randomly or evenly distributed in the solution space. In each 
iteration, particles follow two “extreme values” to update themselves. 
The first value is the optimal solution that the particle found, also 
called the personal best value (pbest); the second extreme value is the 
optimal solution discovered by the entire swarm, also called the global 
best value (gbest). After these optimal solutions are found, the 
particles update their velocity and position based on the following 


equations, as shown in Fig. 18: 
vk+1 = wyt +c, rand} (pbest! — x‘) + cyrand3(gbest* — xs) (16) 


k+1 ky yk+1 
xk+ = xk 4 yk+ (17) 


Where vi is the velocity of particle i of the k™ iteration; w is the 


weight value; c4 and cz are acceleration coefficients (also called 


y 


Fig. 18. Concept of updating the searching path for PSO. 


learning factors), which were used to direct the search toward the 
local and global optimal solutions, respectively; and rand; and 
randz are random numbers between [0, 1]. 

Of all the EC methods, the PSO algorithm is most widely 
applied to GMPPT. Because PSO is an optimization method based 
on swarming, it can conduct GMPPT in distributed architecture 
[54-56] or centralized architecture [57-60]. Regarding algorithms, 
Chowdhury and Saha [54] used adaptive perceptive PSO (APPSO), 
whereas the remaining references used basic PSO algorithms. 
Because implementing APPSO is complex, Ref. [54] provided only 
simulated results, whereas the remaining references provided 
experimental results. This observation shows that the PSO algo- 
rithm is an EC method that can be readily implemented using 
microcontrollers. In terms of centralized architecture, the PSO method 
can accurately track the GMPP in various shading patterns, conduct 
direct cycle control [57-60], and attain reduced steady-state oscillation 
[58]. In [60], random number in the accelerations coefficient is 
removed, developing a deterministic PSO (DPSO) mechanism that 
improved the tracking speed. In addition, Liu et al. [59] proposed 
suggestions regarding practically implementing PSO. 

Conference papers [61-64] involved similar PSO methods; 
authors in [63] and [64] combined PSO and Bayesian fusion, 
rendering a system that can yield excellent GMPPT results in a 
changing environment. 


3.1.3.7. Chaos search method. CS is a stochastic search method 
based on chaos theory, using the ergodicity of the chaos variable 
to conduct an optimization search. Consequently, the search 
results of this method are superior to those of other search 
methods using pure random numbers. Zhou et al. [65] used CS 
method based on dual carrier to conduct GMPP search, employing 
(18) and (19) to randomly produce control variables, measure 
output power, and find the GMPP: 


x! =a+x,(b—a) (18) 


Yi= a+ 4On+2b—@) (19) 


According to [65], CS can accurately find the GMPP under PSC; 
however, simulated and experimental results did not verify this 
theory. 


3.1.3.8. Tabu search method. The Tabu search method used records 
of previous search results to avoid falling into local optimal 
solutions. Zheng et al. [66] proposed using tabu search to solve 
GMPP search problems; however, the authors did not conduct an 
in-depth exploration of its implementation or design guidelines. 
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3.1.4. Other methods 
This subsection introduces additional GMPPT techniques that 
were not categorized into the three previous subsections. 


3.1.4.1. Extreme seeking control method. The extreme seeking 
control (ESC) method is an adaptive close-loop control method 
used to search for extreme values. Compared with traditional P&O, 
ESC enhances tracking speed. Fig. 19 shows the block diagram of a 
typical ESC. To address PSC, Lei et al. [67] used a method similar to that 
of [30] to sweep and search for the GMPP location. The largest 
difference between this method and that of [30] is that it replaced 
P&O by using ESC. Heydari-doostabad et al. [68] improved the tracking 
speed of ESC and the startup power drop, employing a method similar 
to that of [28] to rapidly move the OP close to the GMPP. This method 
exhibits satisfactory tracking speed; however, ESC calculation is 
complex, the internal high-pass filter is susceptible to noise 
interference, and actual implementation is difficult. In addition, ESC 
replaced P&O to provide the control signal required for the power 
switch, hindering ESC integration into traditional PGS firmware. 
Consequently, this method is limited to simulation in current 
literature, and actual experimental data cannot be provided. 


3.1.4.2. Model predictive control method. The model predictive 
control (MPC) method can be used to forecast future events and 
implement controls based on those forecasts; therefore, MPC 
involves optimizing based on the control quantity required in 
the current situation. Bouilouta et al. [69] combined MPC and 
large-step GMPP search to address PSC problems. The MPC 


Fig. 19. The block diagram of ESC. 


Table 1 
Comparison of firmware-based methods. 
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method exhibits favorable tracking speed and accuracy, although 
complex calculations are required. Furthermore, the MPC design 
parameters are related to the system model and the power 
converter topology. Compared with MPC, the design of other 
methods need not take the power converter topology into 
account; thus, actual implementation of MPC is system dependent 
and difficult. Consequently, work in [69] only provided simulated 
results. 


3.1.5. Comparing firmware-based methods 

In this subsection, a simple comparison will be conducted, 
incorporating the journal papers that used the aforementioned 
firmware-based methods suitable for use in centralized architec- 
ture. The comparison addresses the dependence on system archi- 
tecture, tracking speed, probability of successfully finding the 
GMPP, complexity of implementation, and sensors required. 
Table 1 summarizes the comparison results, providing a reference 
for PGS-designers and researchers during method selection. 


3.1.6. Methods for detecting partial shading condition 

The method for detecting PSC is often overlooked in GMPPT- 
related literature mainly because the transient state of various 
operating conditions need not be considered (e.g., operating 
conditions changing from uniform insolation to PSC or vice versa) 
when writing a paper. However, PGS design engineers must be 
able to detect whether the current system is exposed to uniform 
insolation or PSC, and whether the system operating condition 
changed. Because MPPT techniques for handling PSC problems 
must address numerous factors, such techniques are typically 
slower compared with traditional, rapid MPPT methods (e.g., the 
variable-step P&O method). Therefore, an accurate understanding 
of the current operating status facilitates selecting suitable algo- 
rithms and avoiding extending the tracking time. In this subsec- 
tion, the PSC detection methods proposed by the references listed 
in Table 1 are compared. Table 2 summarizes the comparison 
results, showing whether the proposed detection methods can 
accurately distinguish between uniform insolation and PSC, and 
whether they can determine if the GMPP changed and subse- 
quently reinitialize tracking for the five possible operating condi- 
tion changes. Fig. 20 shows the five possible changes in operating 
condition tested in this study. 


Method Dependence Tracking speed Hit count 
on system 
[26-28] Y Fast Low 
[30 N Slow Medium 
[31] N Medium Medium to high 
(relate to alpha) 
[32],[33] N Fast Medium to high 
[36 N Fast High 
[37],[38] N Fast High 
[43], Y Fast Relate to the completeness 
[44] of the training data 
[47] N Fast Medium to high 
[49 N Slow to medium High 
[53 N Slow to medium High 
[57-60] N Slow to medium High 
[67 N Slow to medium High 


Algorithm complexity Sensors required Notes 
Low V,I 
Medium V 
Low V,I Can easily be integrated to original PGS 
Low V,I Can easily be integrated to original PGS 
Medium V,I Requires recursive operation 
Medium V,I Requires recursive operation 
Low to medium Irradiance 
Temperature 
Medium V,I 
Medium to high V,I Stochastic process 
Medium to high V,I Stochastic process 
Medium V,I Stochastic process 
Medium to high V,I Stochastic process 


Table 2 


Comparing PSC detection methods. 
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Parameter description Case | Case II Case III Case IV Case V 


Reinitial Uniform/Non- Reinitial Uniform/Non- Reinitial Uniform/Non- Reinitial Uniform/Non- Reinitial Uniform/Non- 


uniform uniform uniform uniform uniform 
I(Sx4.1): The current of the (k+1)-th Yes No Yes No Yes No Yes No No No 
iteration 
I(S,): The current of the k-th iteration 
AI: A pre-defined constant 
P(k+1): The power of the (k+1)-th Yes No Yes No Yes No Yes No No No 


iteration 


P(k): The power of the k-th iteration 

AP: A pre-defined constant 

V(k+1), P(k+1): The voltage/power of Yes No Yes No Yes No Yes No No No 
the (k+1)-th iteration 


Vk), P(k): The voltage/power of the k-th 

iteration 
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the (k+1)-th iteration 


P(k): The power of the k-th iteration 

AP: Pre-defined constants 

Vailai: The voltage/current sample of Yes Yes Yes Yes Yes Yes Yes Yes No No 
point 1 (in the middle) 


Vazla2: The voltage/current sample of 

point 2 (at the left side) 

Va3,la3: The voltage/current sample of 

point 3 (at the right side) 

d: The searching direction Yes Yes Yes Yes Yes Yes Yes Yes No No 


P,(X2_ 4): New obtained power 
Pk-1(X1 +a): Previous obtained power 
P';(X2_a): Averaged power 

r: A pre-defined constant 


[30,33] AP: The variation of the power Yes No Yes No Yes No Yes No No No 
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3.2. Hardware-based GMPPT techniques 


In addition to using firmware-based methods to address the 
multiple peak values caused by PSCs, many researchers have 
proposed using hardware methods to solve this problem. These 
methods are simply discussed as follows. 


3.2.1. Investigation on system configuration and reconfiguration 
technique 

If unmanaged, PSC cause the useable power of a system to 
decrease. Works in [70-72] studied the three configuration meth- 
ods used in large PGSs, including series-parallel (SP), bridge-linked 
(BL), and total-crossed-tied (TCT). The results showed that in PSC, 
TCT and BL attain superior performance compared with SP. 
Comparing TCT and BL, TCT produces more power. Furthermore, 
works in [73] and [74] proposed immediately adjusting the 
configuration of PV arrays according to the shading pattern to 
compensate for the power losses caused by PSC. Although recon- 
figuration method can mitigate the effects of PSC, this technique 
requires using a switch matrix to implement architecture changes; 
therefore, the system is expensive and the controller design is 
complex. In addition, reconfiguration technique cannot guarantee 
the successful tracking of GMPP in all shading patterns. Ref. [75] 
and Ref. [76] are conference papers discussed system configura- 
tions, whereas Ref. [77-81] are conference papers explored using 
the reconfiguration method to reduce the effects of PSC. 


3.2.2. Multi-level converter 

When a system adopts centralized architecture as shown in 
Fig. 4(a), a multi-level converter can be used when each PV 
module in the PGS must be independently controlled [82-84]. 
Multi-level converter can address PSC problems; however, because 
its number of levels is proportional to the number of series- 
connected PV modules in a PGS, after the architecture of PGS 
changes, the entire hardware system and controller must be re- 
designed. In addition, the greater the number of series-connected 
modules is, the more power semiconductor devices are required, 
increasing costs and power losses. Ref. [85-89] are conference 
papers explored multi-level converters. 


3.2.3. Distributed architecture 

Fig. 4(b)-(d) show that in addition to using a centralized 
architecture, a PGS can be designed to apply a DC/DC converter 
or DC/AC inverter for each PV module. This architecture is called a 
distributed architecture. Because each power circuit has its own 
designated controller, even if PV modules in the system experience 
various levels of shading, the distributed architecture can continue 
to track the GMPP of the entire system. Typically, when a DC/DC 
converter is attached to the back of the PV module as shown in 
Fig. 4(b) and (c), this architecture is called a module-integrated 
converter (MIC) [90-96]. Conversely, when a DC/AC inverter is 
attached to the back of the PV module as shown in Fig. 4(d), this 
architecture is often called a micro-inverter [97]. Distributed 
architecture can guarantee tracking of the system GMPP; further- 
more, because its control is distributed, the reliability and robust- 
ness of the system increase. However, the hardware costs of 
distributed architecture are often high, and implementing the 
controller is complex and requires a high processing speed 
microcontroller. In addition, because micro-inverters are directly 
connected to the utility, the hardware design must address 
grid-interconnection and safety problems [97]. Ref. [98-104] are 
conference papers studied MICs, whereas Ref. [105-107] are 
conference papers explored micro-inverters. 
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Fig. 20. The five possible changes in operating mode tested in this study. 
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Fig. 21. Equalizer-assisted architecture. 


3.2.4, Equalizer-assisted architecture 

Equalizer-assisted architecture as shown in Fig. 21 has been a 
popular method of addressing PSC problems in the past two years. 
Fig. 21 shows that the operating principle of this method is similar 
to the parallel-connected MIC shown in Fig. 4(b). The primary 
difference is that parallel-connected MIC must process all the 
power supplied by the PV module, whereas the equalizer-assisted 
architecture must only process the power difference between 
unbalanced PV modules. Therefore, this architecture guarantees 
tracking of the system GMPP and its hardware costs are less 
compared with that of the parallel-connected MIC method [108- 
112]. Because this method is similar to the equalizer used in the 
series-connected batteries of electric vehicles, researchers have 
called this architecture a PGS equalizer. The advantages and 
disadvantages of this method are the same as those of distributed 
architecture. Conference paper [113] also examined equalizer- 
assisted architecture. 


3.2.5. Summary and comparison of hardware methods 

Previous descriptions show that numerous types of hard- 
ware architecture can solve or mitigate PSC problems. Because 
hardware-based solutions relate to the selection of power con- 
verters, which must be designed to address system specifications 
and power ratings, and most system specifications used in 
literature are distinct. Therefore, no further in-depth discussions 
are conducted regarding hardware-based solutions; however, 
relevant information regarding such solutions is summarized in 
Table 3. 


4. Simulated results of selected methods 


The number of references and the number of citations demon- 
strated that the most commonly used firmware-based techniques 
for addressing PSC were the two-stage, segmental search, and SC 
methods. In this subsection, five representative methods were 
selected and MATLAB was employed to test these methods and 
perform comparisons. The methods selected in this study were 
Ref. [28] (representing the system characteristic curve method, 
cited 73 times), Ref. [30] (the most cited method, cited 215 times), 
Ref. [32] (representing the two-stage searching method, cited 19 
times), Ref. [36] (representing the segmental search method, cited 
55 times), and Ref. [58] (representing the SC method, cited 45 
times). Because the parameter settings of [30] and [32] influenced 
the GMPPT performance, simulations of [30] and [32] were 
conducted using different parameter settings to assess the tracking 
performance. Table 4 lists the various methods and relevant 
parameter setting values. In Section 2, it was demonstrated that 
distinct shading patterns produce distinct P-V characteristic 
curves. However, most references tested only a few specific 
shading patterns; thus, the proposed methods cannot be guaran- 
teed suitable for all possible shading patterns. In this study, 
simulations were performed based on a 5s1p system constituting 
five Sun Power E19/240 PV modules, and MATLAB was used to 
simulate 100-1000 W/m? of irradiance levels on each module in 
100 W/m? increments. Consequently, each module had 10 possible 
irradiance levels. To produce a P-V characteristic curve that 
exhibited five peak values in a 5s1p system, the irradiance level 
for each module must be distinct. In five series-connected modules 
that exhibit different irradiance levels, there are C(10,5)=252 
possible shading patterns. Therefore, the eight methods listed in 
Table 4 were tested using all 252 possible shading patterns in this 
paper. Table 5 lists the results. The specification of the PV module 
used in this study is summarized in Table 6. 
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Table 3 
Comparison of hardware-based methods. 
Method Advantage 


Investigation on system configuration [70-72] 
Reconfiguration methods [73-74] 


No hardware required 
No power converter required 


Multi-level converter [82-84] MPPT for individual PV module 


Distributed power structure (series/parallel MIC) MPPT for individual PV module 
[90-96] High GMPPT accuracy 
System upgradable 


Distributed power structure (micro-inverter) [97] MPPT for individual PV module 
High GMPPT accuracy 
System upgradable 
Compact size 


MPPT for individual PV module 
High GMPPT accuracy 


Balancer auxiliary structure [108-112] 


Only process the power difference between mismatched PV 


modules 
System upgradable 
Compact size 


Table 4 
The methods simulated in this study and their corresponding parameter settings. 


Method Ref. Initial Settings 
Method 1 36 P&O perturbation step in 2nd stage: 1.0 V 
Method 2 30 Search step in 1st stage: 0.6 Vocone 

P&O perturbation step in 2nd stage: 1.0 V 
Method 3 30 Search step in 1st stage: 0.7 Vocone 

P&O perturbation step in 2nd stage: 1.0 V 
Method 4 30 Search step in 1st stage: 0.8 Vocone 

P&O perturbation step in 2nd stage: 1.0 V 
Method 5 32 Search step in 1st stage: 20 W 

P&O perturbation step in 2nd stage: 1.0 V 
Method 6 32 Search step in 1st stage: 50 W 

P&O perturbation step in 2nd stage: 1.0 V 
Method 7 28 NA. 
Method 8 59 o=0.4, cy=1.2, c2=1.6 


In Table 5, hit_count refers to the number of successful times the 
GMPP was located during the 252 tests. Because the second stages of 
[28], [30], and [32] used P&O (locating MPP using P&O yields 
oscillations), tracking success was defined as the located MPP being 
99.8% larger than the actual GMPP; otherwise, the result was regarded 
as tracking failure. In certain shading patterns, the difference between 
the highest and second highest peaks in the P-V characteristic curve 
was small (such as shading pattern [800 W/m?, 400 W/m?, 300 W/m?, 
200 W/m?, 100 W/m], GMPP= 203.9 W, second highest peak= 195.0 
W); this can cause errors in certain GMPPT methods. However, the 
difference in the obtained power will be small if the system tracks the 
second highest peaks instead of the highest ones. Consequently, 
obtaining results that presented a low hit_count and a relatively high 
output power was likely. Therefore, tracking accuracy should also be 
taken into account for performance evaluation. Table 5 lists a 
comparison of tracking accuracy for the eight methods. In Table 5, 
avg_acc refers to the average accuracy value of MPP tracking for the 
252 tests; that is, 252 ,MPPacen/252, and MPPaccn is defined as the 
ratio of the MPP found by the algorithm divided by the actual GMPP in 
the n” test. Similarly, max_acc and min_acc refer to the maximal and 


Disadvantage 


Can only mitigate PSC effect 
Cannot guarantee GMPP tracking 
High system cost 

Complicated control 


Individual PV modules may not operate at their MPP. 
System should be re-designed when PGS architecture 
changes 

High switching loss 


Higher hardware costs 
Complicated control 


Higher hardware costs 

Complicated control 

Grid-interconnection issues should be taken into 
account 


Higher hardware costs 


Complicated control 


minimal values of the 252 test results, respectively. In addition to 
tracking accuracy, tracking speed is also critical to system design. 
When comparing tracking speeds of existing methods, many research- 
ers have compared actual tracking times (unit in seconds); however, 
unless the PGS specifications and power converter of two systems are 
identical (indicating that the hardware settling time is consistent), this 
comparison provides no value. For example, the practical design of an 
MPPT system must address the settling time of power circuits to select 
the interval between two sampling times. If the sampling speed is 
faster compared with the system settling time, the value of the voltage 
and current sampled by the algorithm does not indicate steady state, 
and thus cannot be used to assess the effects generated by the control 
command. However, the settling time of a low-power switching 
converter is shorter than that of a high-power one. Therefore, despite 
using identical GMPPT algorithms, the tracking time of a low-power 
system may be shorter compared with that of a high-power system. In 
other words, assessing the tracking speed of MPPT methods by 
comparing tracking times but without considering other factors 
produces biased results. Therefore, the authors of this paper propose 
that the number of iterations required (related to only the method) 
rather than the time (related to the system and the method) should be 
used to compare tracking speeds. Iteration refers to the number of 
times a specific function is calculated. Using the P&O method as an 
example, one iteration refers to one perturbation. In another example, 
the number of iterations is the number of MPP calculation interrupts 
when a proposed method is implemented by microcontrollers. In 
Table 5, avg_step refers to the average number of required iterations 
for tracking the MPP in the 252 tests. Similarly, max_step and 
min_step refer to the maximal and minimal required iteration 
numbers in the 252 MPPT tests, respectively. Of the eight testing 
methods, only Method 8 (PSO algorithm) exhibited stochastic char- 
acteristics; that is, the tracking results obtained in each execution 
could vary. Thus, test results of Method 8 provided in Table 5 is the 
average of 20 repeated tests using the same shading pattern; conse- 
quently, the max_step and min_step of Method 8 are not integer 
numbers. 
From Table 5, the following can be concluded: 


1) P&O, IncCON, and variable-step P&O/IncCON are the most 
commonly used methods in the second stage of the two- 
stage method. Therefore, the performance of the two-stage 
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Table 5 
Simulated results. 


Direct Patel_0.6 Patel_0.7 Patel_0.8 Power_20 Power_50 Real PSO 
avg_MPPT_acc (%) 98.99 51.36 51.35 98.41 99.85 99.45 52.09 99.96 
max_MPPT_acc (%) 100.00 100.00 100.00 100.00 100.00 100.00 99.79 100.00 
min_MPPT_acc (%) 92.48 17.43 17.43 68.51 96.69 81.21 31.45 99.20 
hit_count 183 11 11 153 227 215 1 251.85 
avg_MPPT_step 20.94 71.44 71.87 473.15 33.98 23.58 42.26 129.19 
max_MPPT_step 31.00 102.00 112.00 1185.00 53.00 38.00 47.00 700.00 
min_MPPT_step 16.00 60.00 65.00 38.00 19.00 13.00 22.00 7.00 


Table6 avoided. The simulated results of this method indicated that 
ae of the sun power E19/240 PV when large steps were selected for the sweep (AP=50 W), the 


Parameters for model 


system required decreased time to find the GMPP (23.58 steps), 
but could miss the GMPP (hit_count=215). Conversely, when 
small steps were selected (AP=20 W), the system had a high 


Maximum power (Pmax) ZOW probability of finding the GMPP (hit_count=227), but required 
Open circuit voltage (Voc) 48.6 V latively | ti 33.98 st Therefi n thi 

MPP voltage (Vmp) 405V a relatively long time ( .98 steps). Therefore, using his 
Short circuit current (Ise) 630A method requires addressing the tradeoff between step size 
MPP current (Imp) 5.93 A and search time. Because this method is simple to implement 


Temperature coefficient (ay) — 132.5 mV/K 


method is determined mainly based on the first-stage (GMPP 
location search). 

2) The probability of finding the GMPP by using a system 
characteristic curve method, which is based on PV panel 
equivalent resistance, is relatively low. In [28], when PSC 
yielded a complex P-V characteristic curve (increasing the 
number of LMPPs), the probability of this method finding the 
GMPP significantly decreased (1/252). Thus, this method 
should only be used for simple P-V characteristic curves (when 
used to assess a 2s1p system, this method can yield a 20/45 
probability of success with avg_acc= 85.42%). 

3) The method used in [30] attained the highest number of citations 
in current GMPPT paper. However, this method is designed after 
observing the P-V characteristic curve caused by PSC and convert- 
ing such observations into simplified search rules. Thus, the 
probability of finding the GMPP is also low. The performance and 
design considerations of this method are outlined as follows: (1) in 
[30], the author recommended using a 0.6-0.7 Vocone step to 
conduct jump search. However, extensively testing by using a 5s1p 
system assembled from Sun Power E19/240 PV modules showed 
that 0.8 Vocone yielded better results. Therefore, prior to using this 
method, simple analysis on system characteristics or parameter 
setting experiments should be conducted to determine the optimal 
design parameters; (2) Although the hit_count (153/252) for this 
method was not high, the avg_acc result (98.41%) indicated a high 
level of tracking accuracy, suggesting that this method did not track 
the real GMPP, but the power of the tracked LMPP was not 
substantially different compared with the power of the GMPP. In 
other words, from the perspective of electricity sales profit, this 
method is worth adopting; and (3) Table 5 shows that large 
numbers of steps are required in this method (473.15), indicating 
that the tracking time is relatively long. Because this method uses 
P&O to track each LMPP, improving this procedure can enhance 
tracking speed. 

4) The method in [32] is another representative method of the 
two-stage searching method. The simulated results showed 
that this method attained excellent hit counts, tracking accu- 
racy, and tracking time. The primary difference between this 
method and that used in [30] was that this method involves 
only a brief scan of the P-V characteristic curve in the first 
stage; therefore, the slow tracking speed encountered in [30] 
(using P&O to track each LMPP reduces tracking speed) can be 


and can be readily integrated into original PGS firmware, it is a 
method worth using, implying two directions of future study. 
First, finding a sweeping rule that can improve hit_count, and 
second, improving the P&O method used in the second stage to 
enhance the GMPPT speed (such as using variable-step P&O). 


5) The method in [36] can also be used to balance hit count, tracking 


accuracy, and tracking time; however, this method cannot be easily 
integrated with original PGS firmware (typically requiring com- 
plete algorithm re-write), making it only suitable for use in newly 
developed PGSs. Although the mathematical principle of this 
method is relatively simple and comprehensive, actual coding 
yields challenges. First, the jump range for sampling points used 
in this method (Va; > Vaz > Va3, See Fig. 13) is the largest of the five 
methods tested, indicating that the firmware designer must 
consider the hardware settling time to avoid excessive stress on 
the power switch. In addition, this method involves using the 
recursive method to gradually reduce the search range; this 
challenges some firmware engineers who are only accustomed to 
sequential programming. 


6) The method in [58] represents the SC method. According to the 


simulated results, the tracking accuracy of this method was 
100% (hit_count=252). Because PSO can be applied to pro- 
blems more complex than those examined in this study (e.g., 
the traveling salesman problem), this result was unsurprising. 
Additionally, the avg_acc of this method was the highest of the 
tested methods, indicating that this method yielded the highest 
electricity sales profit. However, based on the simulated results, 
this method required a relatively long tracking time; therefore 
is unsuitable for use in areas that exhibit rapid shading pattern 
changes. Because PSO is a stochastic optimization method, the 
number of iterations required to track the GMPP for each 
execution may be different; this distinguishes SC from the other 
methods. Furthermore, the parameters in (16) and (17) substan- 
tially affect the tracking performance when this method is used. 
Consequently, optimizing these parameters is a topic worth 
studying. According to [59] and [60], SC can be implemented by 
employing low-cost microcontrollers such as TMS320 or dsPIC33 
series digital signal controllers. 


5. Conclusion 


Extensive recent literature shows that the GMPPT methods 


suitable for use in PSC remain a popular research topic. Current 
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literature comprises numerous algorithms and hardware architec- 
tures related to addressing PSC problems; therefore, PGS designers 
may be challenged to select appropriate firmware or hardware 
architecture. In this study novel GMPPT techniques from diverse 
references were collected, classified, compared, and summarized; 
the advantages and disadvantages of these techniques are dis- 
played in tables. Regarding firmware-based techniques, compar- 
isons were conducted and the indices that must be considered in 
the design were summarized such as dependence on system 
architecture, tracking speed, the probability of successfully track- 
ing the GMPP, algorithm complexity, and required sensors. The 
advantages and disadvantages of the hardware-based techniques 
were also discussed. Furthermore, PSC detection and determina- 
tion methods, which have seldom been discussed in current 
literature, were reviewed and compared. Finally, MATLAB was 
used to extensively simulate five commonly used GMPPT algo- 
rithms and explore their tracking performance; recommendations 
were then proposed for using these methods. Based on the 
simulated results, topics worth studying were also suggested. 
Firmware-based GMPPT methods can continue to be enhanced, 
particularly in terms of maintaining hit count and tracking 
accuracy while accounting for tracking speed, as well as addres- 
sing problems related to uniform insolation and PSC and detecting 
the transient state in an actual system. To conclude, this study 
should serve as a valuable reference for PGS researchers and 
application engineers. 
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